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. A mathematical model for fluid transport in peritoneal dialysis is con- 

structed. The model is based on a three-component nonlinear system of 
two-dimensional partial differential equations for fluid, glucose and albumin 
transport with the relevant boundary and initial conditions. Non-constant 
steady-state solutions of the model are studied. The restrictions on the pa- 
rameters arising in the model are established with the aim to obtain exact 

OO . 

formulae for the non-constant steady-state solutions. As the result, the ex- 
act formulae for the fluid fluxes from blood to tissue and across the tissue 
were constructed together with two linear autonomous ODEs for glucose and 
albumin concentrations. The analytical results were checked for their applica- 
bility for the description of fluid-glucose-albumin transport during peritoneal 
dialysis. 
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1 Introduction 

Peritoneal dialysis is a life saving treatment for chronic patients with end stage 
renal disease (Gokal R and Nolph 1994). Dialysis fluid is infused into the peri- 
toneal cavity, and, during its dwell there, small metabolites (urea, creatinine) and 
other uremic toxins diffuse from blood to the fluid, and after some time (usually 
a few hours) are removed together with the drained fluid. The treatment is re- 
peated continuously. The peritoneal transport occurs between dialysis fluid in the 



1 



peritoneal cavity and blood passing down capillaries in tissue surrounding the peri- 
toneal cavity. The capillaries are distributed within the tissue at different distance 
from the tissue surface that is in contact with dialysis fluid. The solutes, which are 
transported between blood and dialysis fluid, have to cross two transport barriers: 
the capillary wall and a tissue layer. Typically, many solutes are transported from 
blood to dialysate, but some solutes that are present in high concentration in dial- 
ysis fluid are transported to blood. This kind of transport system happens also in 
other medical treatments, as local delivery of anticancer medications, and some ex- 
perimental or natural physiological phenomena. Mathematical description of these 
systems was obtained using partial differential equations based on the simplifica- 
tion that capillaries are homogeneously distributed within the tissue (Flessner et al. 
1984; Waniewski et al. 1999; Waniewski 2002). Experimental evidence confirmed 
the good applicability of such models (Flessner et al. 1985). 

Another objective of peritoneal dialysis is to remove excess water from the patient 
(Gokal R and Nolph 1994). This is gained by inducing high osmotic pressure in 
dialysis fluid by adding a solute in high concentration. The most often used solute 
is glucose. This medical application of high osmotic pressure is rather unique for 
peritoneal dialysis. Mathematical description of fluid and solute transport between 
blood and dialysis fluid in the peritoneal cavity has not been formulated fully yet, in 
spite of the well known basic physical laws for such transport. A previous attempt 
did not result in a satisfactory description, and was disproved later on (Seames et 
al. 1990;Flessner 1994). Recent mathematical, theoretical and numerical studies 
introduced new concepts on peritoneal transport and yielded better results for the 
transport of fluid and osmotic agent (Cherniha and Waniewski 2005; Stachowska- 
Pietka et al. 2005; Stachowska-Pietka et al. 2006; Cherniha et al. 2007; Waniewski et 
al.2007, Waniewski et al.2009). However, the problem of a combined description of 
osmotic ultrafiltration to the peritoneal cavity, absorption of osmotic agent from the 
peritoneal cavity and leak of macromolecules (proteins, e.g., albumin) from blood to 
the peritoneal cavity has not been addressed yet, see for example in (Flessner 2001; 
Stachowska-Pietka et al. 2007). Therefore, we present here a new extended model 
for these phenomena and investigate its mathematical structure. 

The paper is organized as follows. In section 2, a mathematical model of glu- 
cose and albumin transport in peritoneal dialysis is constructed. In section 3, non- 
constant steady-state solutions of the model are constructed and their properties 
are investigated. Moreover, these solutions are tested for the real parameters that 
represent clinical treatments of peritoneal dialysis. The results are compared with 
those derived by numerical simulations for a simplified model (Cherniha et al. 2007). 
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Finally, we present some conclusions and discussion in the last section. 

2. Mathematical model 

The mathematical description of transport processes within the tissue consists in 
local balance of fluid volume and solute mass. For incompressible fluid, the change 
of volume may occur due to elasticity of the tissue. The fractional void volume, 
i.e. the volume occupied by the fluid in the interstitium (the rest of the tissue being 
cells and macromolecules forming interstitium) expressed per one unit volume of the 
whole tissue is denoted v(t,x), and its time evolution is described by the following 
equation: 

dv dju 

where ju(t,x) is the volumetric fluid flux across the tissue (ultrafiltration), qu(t,x) 
is the density of volumetric fluid flux from blood to the tissue, while the density 
of volumetric fluid flux from the tissue to the lymphatic vessels q\ (hereafter we 
assume that it is a known positive constant, nevertheless it can be also a function 
of hydrostatic pressure) produces absorbtion of solutes from the tissue. 

The independent variables t is time, and x is the distance from the tissue surface 
in contact with dialysis fluid (flat geometry of the tissue is here assumed). The 
solutes, glucose and albumin, are distributed only within the interstitial fluid, and 
their concentrations in this fluid are denoted by Ca(t,x) and Ca(£, x), respectively. 
The equation that describes the local changes of glucose amount, vCq, is as follows: 

where ja(t, x) is glucose flux through the tissue, and qc(t, x) is the density of glucose 
flux from blood. 

Similarly, equation that describes the local changes of albumin amount, uC Al is 
as follows: 

d(avC A ) dj A , /q\ 

where jA(t,x) is albumin flux through the tissue, and qA(t, x) is the density of 
albumin flux from blood, while the coefficient a < 1 takes into account that a part 
of the fractional void volume v is only accessible for albumin because of big size of 
molecules (Flessner 2001; Stachowska-Pietka et al. 2007). The flows of fluid and 
solutes are described according to linear non-equilibrium thermodynamics. Osmotic 
pressure of glucose and oncotic pressure of albumin are described by van't Hoff law, 
i.e. it is proportional to the relevant concentrations. 
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The fluid flux across the tissue is generated by hydrostatic, osmotic and oncotic 
pressure gradients: 

ju = -vK— + a TG vKRT— ° + o TAl vKRT—^, (4) 
ox ox ox 

whereas for the density of fluid flux from blood to tissue we assume that it is gen- 
erated by the hydrostatic, osmotic and oncotic pressure differences between blood 
and tissue: 

qu = L p a(P B - P) - L p aa G RT(C GB - C G ) - 1 L p aa A RT(C AB - C A ), (5) 

where P(t, x) is hydrostatic pressure. 

The glucose solute flux across the tissue is composed of diffusive component 
(proportional to glucose concentration gradient) and convective component (pro- 
portional to glucose concentration and fluid flux): 

dC G 

jo = -vD G -^— + S TG C G ju. (6) 

The density of glucose flux from blood to tissue has diffusive component (propor- 
tional to the difference of glucose concentration in blood, C GB , and glucose concen- 
tration in tissue, C G ), convective component (proportional to the density of fluid 
flow from blood to tissue, qu) and component presenting lymphatic absorbtion: 

q G = p G a(C GB - C G ) + S G q v ((1 - F G )C GB + F G C G ) - q x C G . (7) 

In quite similar way, we can construct equations for the albumin solute flux across 
the tissue,j A (t,x), and the density of albumin flux from blood to tissue q A (t,x): 

j A = -auD A ^^- + S TA C A j v , (8) 



q A = p A a(C AB - C A ) + S Aqu ((1 - F A )C AB + F A C A ) - q x C A . (9) 

The coefficients in the above equations are: K - hydraulic permeability of tissue, 
otg an d <7ta ~ the Staverman reflection coefficients for glucose and albumin in tissue, 
Stg = 1 — °tg and Sta = 1 — ctta- sieving coefficients of glucose and albumin in 
tissue, a G and a A - the Staverman reflection coefficients for glucose and albumin in 
the capillary wall, S G = 1 — o G and S A — 1 — o A - sieving coefficients of glucose and 
albumin in the capillary wall, R - gas constant, T - temperature, L p - hydraulic 
permeability of the capillary wall, a - density of capillary surface area, D G and D A 
- diffusivities of glucose and albumin in tissue, p G and p A - diffusive permeabilities 
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of the capillary wall for glucose and albumin, Pg - hydrostatic pressure in blood, 
Fq < 1 and Fa < 1 - weighing factors, and 7 < 1 is a coefficient that recalculates 
osmotic pressure of albumin to the total oncotic pressure exerted by all proteins 
(Stachowska-Pietka et al. 2007). 

Equations (Cp)-© together with equations (J4j)— (J9j) for flows form a system of 
three nonlinear partial differential equations with three variables: u, P, Ca, and C G . 
Therefore, an additional, constitutive, equation is necessary, and this is the equation 
describing how fractional fluid volume, z/, depends on interstitial pressure, P. This 
dependence can be established using data from experimental studies ( Stachowska- 
Pietka et al. 2006). It turns out that 

v = F(P), (10) 

where F is a monotonically non-decreasing bounded function with the limits: F — > 
z/ m ; n if P — > P m in and F — > z/ max if P — > P max (particularly, one may take P min = 
— 00, P max = 00). Here v m \ n < 1 and v max < 1 are empirically measured constants. 
The reader may find possible analytical forms for the function F in (Stachowska- 
Pietka et al. 2006; Cherniha et al. 2007). 

Boundary conditions for a tissue layer of width L impermeable at x = L and in 
contact with dialysis fluid at x = follows: 

x = 0: P = P D , C G = C GD , C a = Cad (11) 



,= L:£ = 0, f^O, f^0. (12) 
OX ox ox 

Initial conditions describe equilibrium within the tissue without any contact with 
dialysis fluid at x = 0: 

t = : P = P , C G = C GB , C A = C AB . (13) 

It is easily seen that equations ffT])- f|T0|) can be united into three nonlinear partial 
differential equations (PDEs) for finding the hydrostatic pressure P(t, x), the glucose 
concentration C G and the albumin concentration Ca- Thus, these PDEs together 
with boundary and initial conditions (JTTT)- (fl3l) form the nonlinear boundary- value 
problem. Possible values of the parameters arising in this problem can be established 
from experimental data published in a wide range of papers. 

Finally, we note that finding the function ju(t,x), presenting the fluid fluid flux 
across the tissue and calculating ultrafiltration from the tissue to peritoneal cavity, 
is the most important. 
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3. Steady-state solutions of the model and their applications 

First of all we consider the special case: if the boundary conditions ffTTj) is 
replaced by zero Neumann conditions at x = then the steady state solution can 
be easily found because it does't depend on x. In fact solving algebraic equations 

qu-qi = 0, q G = 0, q A = 0, (14) 

one easily obtains the constant steady-state solution 

_ VgO-+QiSg{1-F g ) (~i 
U G - PG a+ qi {l-SGFG)^ GB 

°A - PA a+ qi (l-S A F A ) L 'AB (15) 



P *= P B-*{-ta+ RT ( P Ga +qi tsG F G) + 



PGO,+qi{l-S G F G ) p A a+qi(l-S A F A ) 

In the case qi = 0, i.e., zero flux from the tissue to the lymphatic vessels, formulae 
( fTBT) produce 



Cq — Cgb, C\ — Cab, P* — Pb, (16) 
otherwise 

C* G <C GB , C* A <C AB , P*<P B - (17) 

However, such solution cannot describe any peritoneal transport. 

To find non-constant steady-state solutions, we reduce Eqs. (CO)-© to an equiv- 
alent form by introducing non-dimensional independent and dependent variables 

x ,^_ (18) 

p(f, **) = u(f, x>) = ° c - C °" , w(t\ x') = C ' A . (19) 

"D — "0 °GD _ ^GB <^GD ~ ^GB 

Thus, after rather simple calculations and taking into account Eq. (BJ, (jBJ), and (JHJ), 
one obtains Eq. ([I])-© in the forms (hereafter upper index * is omitted) 

1 dv Id /vdp\ d fi>du\ d fvdw\ ^ ,^q\ 
to dt t dx V dx J dx\dx) dx\dx) ' 

1 d(vu) ^ aTGUodv ^ d rvdu\ ^ Stg 9 /uudp\ ^ d /uvdu^ 
t dt t dt dx\dx) t dx\dxJ dx\dx) 

-Stg°2-^ ( UU qJ J ) + ( u o( S g - Stg) + S G F G u)qu - b x u - a TG u Q qi, (21) 
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' ^0 


dx 


\ dx ) 



i'du\ 



to dt to dt dx V dx ) t dx V dx ) dx\ dx 

d (w*vdw\ 

-S TA o 2 —y dx J + {w (S A - S TA ) + S A F A w*)qu - b 2 w* - a TA w q h (22) 

where 

& = t (£fa> - p) + ^ + ^*) > 

a x = <j tg KRT Cg °- 2 Ggb , a 2 = <r TA KKT'y CaD ^ B , 

dl = % d 2 = ^ (23) 
&i = Pea + ft, & 2 = PaO + qi, 

to = K(.p L D -p y w*=w-w , 

We want to find the steady state solutions of Eqs. (!20|) - (|22|) satisfying the bound- 
ary conditions f[TTj) -( TT2"j) . They take the form 

x = : p = 1, u = 1, w= ° AD (24) 

L-GD — <^GB 

dp du dw , , 

x = l: / = 0, — = 0, — = 0. 25 
dx dx dx 

for the non-dimensional variables. 

Obviously, Eqs. ( l20|) - (!22!) can be reduced to the ordinary differential equations 
(ODEs) to find steady state solutions: 

1 d (vdr>\ d (vdu\ d (vdw\ , , 

T T* (-£) " ) " hr) + * " 9 = °- < 26 > 



^ d rvdu\ ^ Stg d /uudp\ ^ d /uudu\ 
dx\dxJ t dx\dxJ dx\ dx J 

-S TG a 2 -^-( UP f W S ] + (u (S G - S TG ) + S G F G u)qjj - b x u - (T TG u qi = 0, (27) 
dx \ dx J 

, d fudw\ St A d fw*i>dp\ d fw*i>du\ 

«2-H —r~ + ~, H — ; — - — 1 — - — 



' dx\ dx J t dx V dx ) TA 1 dx V rfx / 
' StaC72 7^{ W ^ W ) +(^o(5'a- ^ta) + S A F A w*)q v -b 2 w* - a TA w qi = 0. (28) 



7 



Non-linear system of ODEs (I26l) - (l28j) is still very complicated for integrations with 
arbitrary coefficients. Thus, we look for the correctly-specified coefficients when this 
system can be simplified. It can be noted that the relations 

Sa = Sta, Sq = Stg (29) 

lead to an essential simplification of this system. In fact, using equation (J26l) . ex- 
pressions for qu from ( |23i) and ju from (jl]), rewritten in non-dimensional variables 

r ( 1 dp du dw\ . . 

JU = Lv{--f- + <T X — + £7 2 —), 30 

V t ox Ox Ox / 
one arrives at the semicoupled system of ODEs 

K d f vdqu \ =qu _ qi ( 31 ) 
L p aL 2 dx\ dx J 

ju = KV d<1U (32) 
LpdL dec 

to find functions qu and ju- 

Equation ( 13TT) is the linear second-order ODE provided the function v is known. 
Nevertheless v depends on the pressure, which is also to be found function, we may 
use additional restrictions on the given function F from fflQ|) . For example, assuming 
that F is proportional to the inverse function to P(x), we obtain 

v{x) = u m , (33) 

where v m is a positive constant. Substituting (|33l) into system (|3TT) -(j32l). we easily 
find its general solution: 

q u = C 1 e- Xx + C 2 e Xx + q l , (34) 



^(-C^ + Cvn X — ]J (35) 

The arbitrary constants C\ and C 2 can be specified using the boundary conditions 
( I24|) -( |25|) since the functions qu and ju are expressed via p, u, w and their first- 
order derivatives (see formulae ( l23l) and ( l30|) ). Making rather simple calculations, 
one obtains 

e 2A 1 

c ^ = ^w)Y^> c '2 = (go-*) T ^ 2T , (36) 
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where 



L n aL 2 ( \ , _ s Cad — C 



a FHffo- 1 )+" + ' V x ° r )• (37) 

Having the explicit formulae for q u and jVj system of ODEs fl27|) and fl28|) with 



restrictions (1291 can be reduced to two linear autonomous ODEs 

d 2 u Sn,~ _v,. „ ^„,du 



and 



^+^(C' 1 e-^-C' 2 e^)^+(/ 2 (C' 1 e-^+C' 2 e^)-/ t2 )( W - W o)-^oi =0 (39) 
dor A dx \ / 

to find the functions and u'(x). Hereafter the notations 

/i = S G F G - 5 G , «! = p G a + (1 - S G F G )q h u 01 = a G u q h 
/ 2 = S'a-Fa - SU, «2 = P^a + (1 - S A F A )qi, w i = (T A w qi 

are used. To our best knowledge, the general solutions of ODEs fl38|) and (139]) (obvi- 
ously, both equations have the same structure) are unknown in explicit form. On the 
other hand, since the unknown functions u(x) and w(x) should satisfy the boundary 
conditions (|24p - (|25p . the corresponding linear problems can be numerically solved 
using, for example, Maple program package. Finally, using two expressions for q v 
from ( 1231 and (134|) . we obtain the function 

p(x) =p + t a lU + t a 2 (w - w ) - (de^* + C 2 e Xx - q^j . (41) 

In the next section, the realistic values of parameters arising in the formulae 
derived above will be established and used to calculate the steady-state solutions of 
the model with restrictions f[2"9l and 



Since restriction ( 133]) is rather an artificial simplification for modeling fluid trans- 
port in peritoneal dialysis and it needs additional justifications, we examined the 
cases when the function v is non-constant and satisfies the properties presented after 
formula ([TO]) . The simplest case occurs when v is linear monotonically decreasing 
function 

v(p{x)) = u(x) = v max - {v max - v min )x, x G [0, 1]. (42) 
Substituting (j4"2l into (1291 . we obtain the linear ODE 

i t \ \d 2( lu , s.dqu L p aL 2 

yVmax — \y m ax ~ ^min)X)-j^ {^rnax ~ v min)—^ {Qu ~ Ql) ~ U - (4oJ 
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It turns out that this ODE reduces to the modified Bessel equation of the zero order 
(see, e.g., Bateman 1974) 

2 d 2 q v dq v 2 
by the substitution 

y = 4:d*(is*-x), q v = qu~Qi, v* = > 1, <)* = —. > (45) 

Umax ^min •*»■ y^max L'min) 

The general solution of (144"]) is well-known, hence, using formulae (115]) . one obtains 
the solution of ( |43l) : 

gc/ = Ci/ (2 V / 5*(^ - x)) + C 2j Rr (2v / 5*(^ - 5) + fl, (46) 

where Jo and K are the modified Bessel functions of the first and third kind, re- 
spectively. 

Substituting the function qu obtained into (I3TT) and using the well-known rela- 
tions between the Bessel functions ( Bateman 1974), we find the function: 



-L^^-^fcj^y/SM-x)) - C 2 K 1 (2 v /^(z/,-x))), (47) 



where I\ and K\ are the modified Bessel functions of the first order. The arbitrary 
constants C\ and C 2 can be specified in the quite similar way as this was done in the 
case v{x) = const. Omitting rather simple calculations, we present only the result: 



c = (90-^(2^(^-1)) 

1 I (2V^)#i(2y^ - 1)) + AT (2v^)/i(2^*(^ - 1))' 

C 2 = ( 9 o-^)/ 1 (2 v ^W)) 

J (2v / &)^i(2 v /^(^ - 1)) + K (2VS^)h(2^/S^ - 1)) 



where go is defined by (J3j 

Thus, we have found the explicit formulae for qu and ju- In contrast to the 
case u(x) = const, they contain not only elementary functions but transcendental 
functions as well. Having formulae (jl6"]) - ()47p . system of ODEs (1271) and (128]) with 
restrictions (129]) can be reduced to two linear autonomous ODEs to find the functions 
u(x) and w(x). These equations possess the forms 

-^±(v^T r ^(Ci/i(2 v /5 + (^ -x)) - C 2 K 1 (2 v /^(z/,-x))) M ) (50) 
+ (s G F G (C 1 I (2 y /5*(v*-x)) + C 2 K (2y/5*(v*-x))) -k^u- u i = 

10 



and 





■s/v* - x(Ci/i(2a/(5*(^* - xj) - C 2 Ki(2y/5*(v* - x)))(w - w )j 




Nevertheless both equations are linear second order ODEs with the same struc- 
ture, it seems to be unrealistic to construct their general solutions because of their 
awkwardness. Thus, we will numerically solve them together with the boundary 
conditions (I2~4"j) -f l2"5"j) . using Maple program package. In the next section, the real- 
istic values of the parameters arising in formulae (H6l)- (l5TT) will be established and 
used to calculate the steady-state solutions. 



Here we present the results of application of the formulae derived in Section 3. 
Our aim is to check whether they are applicable for describing the fluid-glucose- 
albumin transport in peritoneal dialysis. The parameters arising in the formulae 
were derived from experimental data and applied in previous mathematical studies 
(Imhofz et al. 1998; Zakaria et al. 1999; Flessner 2001; Waniewski 2001; Waniewski 
2004; Stachowska-Pietka et al. 2006; Cherniha et al. 2007). 

Thus, we used the following values of parameters and absolute constants: 
K = 5.14-10 -5 cm 2 -mirT x -mmH g~ l - hydraulic permeability of tissue, RT = 18-10 3 
mmHg ■ mmol~ l ■ mL - gas constant times temperature, L = 1.0 cm - width of the 
tissue layer, Lpa = 7.3 • 10~ 5 mL ■ min~ l ■ mmHg" 1 ■ g" 1 - hydraulic permeability of 
the capillary wall, qi = 0.26 • 10 -4 mL ■ min' 1 ■ cm" 3 - volumetric fluid flux from the 
tissue to the lymphatic vessels, Dq = 12.11 • 10 -5 cm 2 -mirC x - diffusivity of glucose 
in tissue divided by z/ , Da = 0.2 ■ 10 -5 cm 2 ■ mirT 1 - diffusivity of albumin in tissue 
divided by i> , pea = 3.4- 10~ 2 mL-min" 1 -g" 1 - diffusive permeability of the capillary 
wall for glucose, pacl = 3-10 -4 mL-min~ x -g~ x - diffusive permeability of the capillary 
wall for albumin (there are no precise values for this parameter in literature, however, 
Pa is expected to be at least 100 times smaller than pc [13] ), otg - the Staverman 
reflection coefficient for glucose in tissue (varies from to 0.01), Stg — 1 — °"tg - 
sieving coefficient of glucose in tissue, gta - the Staverman reflection coefficient for 
albumin in tissue (varies from 0.05 to 0.5), St a = 1 — <Jta - sieving coefficient of 
albumin in tissue, Cqb = 6 ■ 10~ 3 mmol ■ mL" 1 and Cab = 0.6 • 10~ 3 mmol ■ mL" 1 - 
glucose and albumin concentration in the blood, respectively, Cqd = 180 • 10~ 3 
mmol ■ mL -1 and Cad = - glucose and albumin concentration in the dialysate, 



4. Applications for the peritoneal transport 
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0,2 0,4 0,6 0,8 1 0.2 0,4 0,6 0,8 1 



Figure 1: The fluid flux from blood to tissue qu (in mL ■ min 1 • cm 3 ) and the fluid 
flux across tissue ju (in rnL ■ min' 1 ■ cm~ 2 ) in the case v = (v max + u min ) /2 = 0.26, 
c?tg = 0.001, and o TA = 0.25. 

respectively, Pg = 15 mmHg - hydrostatic pressure in the blood, Pp = 12 mmHg - 
intraperitoneal hydrostatic pressure, Pq = - initial interstitial hydrostatic pressure, 
F G = 0.5 and Fa = 0.5 - weighing factors for glucose and albumin, respectively, and 
the non-dimensional coefficients u min = 0.17, v max = 0.35, u = 0.17, a = 0.8,7 — 1- 

Let us consider the first case, when restrictions fl29|) and fl33|) take place. First 
of all, it seems to be reasonable to set v m = [y max + ^ m j„)/2 = 0.26, i.e., we assume 
that the fractional void volume at the steady state stage of the the peritoneal trans- 
port is an intermediate value between its maximum and minimum. The Staverman 
reflection coefficients for glucose and albumin in tissue are now otg — 0.001 and 
Ota = 0.25, respectively. 

Fig.l presents the space distributions of the density of fluid flux from blood to 
tissue qu and the fluid flux across tissue ju, calculated using formulae (13^1) -(1371). One 
sees that the function quix) is monotonically decreasing, while the function ju(x) 
is monotonically increasing with the distance from the peritoneal surface, and this 
corresponds to the experimental data and numerical simulations for the simplified 
model (Cherniha et al. 2007; Waniewski et al. 2007). 

Using the value of the fluid flux ju at the point x = 0, one may calculate the 
reverse water flow (i.e. out of the tissue to the cavity). Total fluid outflow from the 
tissue to the cavity, calculated assuming that the surface area of the contact between 
dialysis fluid and peritoneum is equal to 0.5 m 2 , is about 0.50 mL/min. Note the 
similar value was obtained in (Cherniha et al. 2007) using numerical simulations for 
the simplified model. 
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Figure 2: The non-dimensional glucose concentration u and albumin concentration w 
in the case v = {u max + u min )/2 = 0.26, otg — 0.001, Ota = 0.25, and p^a = 3 • 10~ 4 
mL ■ mirT 1 ■ g^ 1 . 

Fig. 2 presents the space distributions of the glucose and albumin concentrations. 
The glucose concentration u decreases rapidly with the distance from the peritoneal 
cavity to zero in the deeper tissue layer and is practically zero for x > 0.15. As 
one may expect, the albumin concentration w increases with the distance from 
the peritoneal cavity. However, the concentration increases rapidly only in a very 
thin layer, while for any x > 0.05 it is practically constant, w ~ 0.0034, and in 
agreement with the initial albumin profile Ca = Cab = 0.6 • 10~ 3 mmol ■ mL -1 (see 
(j!3p ). However, such sharp increase in concentration is rather unrealistic and is a 
consequence of the assumption (|33|) . 

Let us consider now the second case, which is more realistic. We assume that 
restrictions (1291) and (142p take place. The Staverman reflection coefficients were 
chosen the same as in the first case, i.e., otg = 0.001 and ax a = 0.25. The results 
are presented on Fig. 3 and 4. It is quite interesting that the profiles for the functions 
qu(x) and ju(x) pictured on Fig. 3 are very similar to those on Fig. 1, although the 
relevant formulae are essentially different (the reader may compare ( 146]) - ( 1491) with 
(1341 -( 1371) ). Moreover, the total fluid outflow from the tissue to the cavity, calculated 
under the above mentioned assumption is approximately equal to —0.55 mL/min, 
that is only about 10 percents higher than in the previous case. 

The space distributions of the glucose and albumin concentrations are pictured 
on Fig.4. The glucose concentration u is again a decreasing function, however the 
tissue layer with non-vanishing u is wider. 

The main difference occurs in the case of the albumin concentration w. We 
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Figure 3: The fluid flux from blood to tissue qu (in mL ■ min 1 • cm 3 ) and the fluid 
flux across tissue jv (in rnL ■ min' 1 ■ cm~ 2 ) in the case v = v max — [v max — v min )x, 
a TG = 0.001, and a TA = 0.25. 

found that the albumin concentration w essentially depends on the parameter paQ- 
Three curves pictured on Fig. 4 correspond to the values of the diffusive permeability 
Pao, = 5 • 10~ 4 , 3 • 10~ 4 , and 2 • 10~ 4 mL ■ min^ 1 ■ g" 1 , respectively. Negative values 
occurring close to the peritoneal cavity (the red curve corresponding to the smallest 
Paci value) are very small and can be either a consequence of the simplifications 
in the model or an error in numerical solving ODE ( 15T|) . We may interpret these 
negative values as follows: there is no albumin in this layer of the tissue because it 
was already removed to the peritoneal cavity. 

Finally, one observes that the relevant curves on Fig. 2 and Fig. 4 (the middle 
curve) are essentially different, nevertheless they were obtained for the same param- 
eters. As follows from the Fig.4, the albumin concentration w increases in deeper 
layers of tissue as well, and is equal to the constant, corresponding to the initial 
profile, only close to the opposite side of tissue. Thus, such steady-state profile of 
the albumin concentration is more realistic than in the first case (see Fig.2). 

5. Conclusions 

In this paper, a new mathematical model for fluid transport in peritoneal dial- 
ysis was constructed. The model is based on a three-component nonlinear system 
of two-dimensional partial differential equations and the relevant boundary and ini- 
tial conditions. To analyze the non-constant steady-state solutions, the model was 
reduced to the non-dimensional form. Having in mind to obtain exact formulae for 
such solutions, we found the restrictions on parameters arising in the model that 
essentially simplified the equations of the model. As the result, the exact formulae 
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Figure 4: The non-dimensional glucose concentration u and albumin concentration 
w in the case v = v max -{v max -v min )x, a TG = 0.001, o TA = 0.25, and p A a = 5-10" 4 , 
3 • 10~ 4 , 2 • 10" 4 mL ■ min' 1 ■ g' 1 

for the density of fluid flux from blood to tissue and the fluid flux across the tis- 
sue were constructed together with two linear autonomous ODEs for glucose and 
albumin concentrations. 

The analytical results were checked, whether they are applicable for describing 
the fluid-glucose-albumin transport in peritoneal dialysis. Thus, the realistic values 
of the parameters, arising in the formulae, were used to calculate the steady-state so- 
lutions of the model. The conclusion is rather optimistic because even the simplest 
approximation of fractional fluid volume v via linear function with the correctly 
specified coefficients leads to plausible results. Other, more realistic approximation 
of v may in future result in similar exact formulae. However, in general the assump- 
tion about equality of the reflection coefficients in the tissue and in the capillary 
wall, although it demonstrates an interesting specific symmetry in the equations, 
can be too restrictive for practical applications of the derived formulae ( Waniewski 
et al. 2009). Therefore, other approaches to find the analytical solutions of the 
model need to be looked for. 
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